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Abstract 

Water activity and its relationship with interactions stabilising lamel¬ 
lar stacks of mixed lipid bilayers in their fluid state are investigated by 
means of osmotic pressure measurements coupled with small-angle x-ray 
scattering. The (electrically-neutral) bilayers are composed of a mixture 
in various proportions of lecithin, a zwitterionic phospholipid, and Simul- 
sol, a non-ionic cosurfactant with an ethoxylated polar head. For highly 
dehydrated samples the osmotic pressure profile always exhibits the “clas¬ 
sical” exponential decay as hydration increases but, depending on Simulsol 
to lecithin ratio, it becomes either of the “bound” or “unbound” types for 
more water-swollen systems. A simple thermodynamic model is used for 
interpreting the results without resorting to the celebrated but elusive 
“hydration forces”. 
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1 Introduction 

From mixtures of lipids and water, multilamellar systems may naturally emerge. 
In such self-assembled systems, a periodic structure is formed by stacking lipid 
bilayers and layers of water. In the case where there is no in-plane order in the 
bilayers-in their so-called “fluid state”-, the system exhibits the symmetry of a 
smectic A phase, commonly referred to as a lamellar phase in the context of 
lipid materials [T]. 

Assessing the mechanisms responsible for the stability of such lamellar struc¬ 
tures has been the core motivation for quite a large number of experimental or 
theoretical studies, with “hydration interactions” [2], on one hand, and “undu¬ 
lation interactions” [3] on the other hand emerging as central concepts more 
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than 35 years ago when direct electrostatic interactions are irrelevant. Tech¬ 
niques of choice for studying inter-bilayer interactions include osmotic pressure 
control 111I1I3, line-shape analysis in high-resolution [S] or grazing-incidence 
small-angle scattering geometries [7], dynamic light scattering [S], and “direct” 
methods with, e.g. Surface Force Apparatuses [H]. Clear experimental evi¬ 
dences have been obtained as regards undulation interactions in the swollen 
end of dilution lines [5], while the other (highly dehydrated) limit has often 
been generically characterised in terms of “hydration interactions” from the 
observed exponentially-decaying force (or pressure) profile with natural scale 
0.2-0.4 nm 0. 

Here, we study a lecithin-based lamellar system in the presence of a non¬ 
ionic co-surfactant, a system which has shown its ability to efficiently encapsu¬ 
late DNA fragments in spite of the absence of any obvious direct electrostatic 
mechanism at play nuHH]. The co-surfactant we use is an ethoxylated fatty 
acid, i.e. a (short) non-ionic block copolymer with amphiphilic properties. The 
present work therefore somehow expands Ref. [umi, the focus being now os¬ 
motic pressure control instead of line-shape analysis in small-angle scattering. 

The lamellar structure of the stacked bilayers is equilibrated with various 
semi-dilute (aqueous) solutions of polymers, which gives a handle on water 
activity (or, equivalently, osmotic pressure H), while the stacking period i is 
determined by means of small-angle X-ray scattering, following the method 
popularised by V.A. Parsegian [^. Other relevant quantities are commonly ma¬ 
nipulated when presenting or discussing the results, namely the bilayer volume 
fraction cj) and the interfacial area per (average) amphiphilic molecule S. As¬ 
suming homogeneous and ideally flat bilayers, a simple geometric description of 
the lecithin-Simulsol lamellar stack in water gives (see, e.g., Ref. M) 


e = 2lx- 


( 1 ) 


where v is the effective molecular volume xvs + (1 — x)vl of the bilayer species 
(with vl, respectively vs, being the molecular volume of the lecithin, resp. 
Simulsol molecules and x the Simulsol mole fraction in the bilayer), E similarly 
being an effective interfacial molecular area xY,s -I- (1 — x)J1l derived from ac¬ 
tual Simulsol Es and lecithin E^, interfacial areas. Note that, while H, £, x or 
(j) are experimentally well-defined quantities, and vl, as well as vs, can safely 
be considered as constant parameters, E is model-dependent, as should be clear 
from the assumptions leading to eq. H]) : These assumptions would be spoilt by 
the presence of numerous structural defects (holes across single bilayers, pas¬ 
sages connecting adjacent bilayers-implying no longer homogeneous bilayers), 
or large-amplitude area-storing bilayer undulations-implying no longer flat bi¬ 
layers. 
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2 Thermodynamic considerations 


2.1 Lamellar stacks 

Implicitly assuming for the sake of simplicity a two-component (lipid-water) 
system, it is customary to cast the interpretation framework of the n(£) data 
into the mould of, loosely speaking, an inter-bilayer interaction potential energy 
per unit bilayer area V, but perhaps more rigorous (as already noticed in Ref. 0) 
to start from the excess free energy of the bilayer stack. Per unit volume of the 
lamellar stack (and disregarding an obvious dependence on temperature that 
remains implicit in the following), the excess free energy is a function of two 
among the three quantities appearing in eq. ([T]). Choosing (jj and E, the excess 
free energy density, formally written as follows: 

/exc(<(>,S) (2) 


should be zero, by definition, for either </> = 0 or ^ = 1. The total Gibbs free 
energy is then expressed as 


G — {NjjjVw Niv) 


p _(_ ^^(1 _ ^) _|_ ^^(j) + /exc(<(>, S) 
Vw V 


(3) 


Here, p is the (hydrostatic) pressure. Indices w, respectively I, refer to water, 
resp. “average” lecithin / Simulsol lipid. The total numbers of molecules of a 
given species are N^i and Ni , and the chemical potentials for the pure species are 
and . The lipid volume fraction \s (p = Niv/{N^Vw + Niv). Incompress¬ 
ibility is enforced by considering the molecular volumes Vw and v as constant 
parameters, with therefore dG/dp = N^Vw + Niv. 

The “dilution law”, namely the ^{(p) relation, is derived from eq. ([T]) once 
the optimal interfacial area is obtained, for a given composition cp, by solving 
the minimisation equation 

which is equivalent to solving 


The “bilayer equation of state”, that is to say the explicit solution to this 
latter equation, is commonly expressed as E((/)) or E(I') relations. 

The “lamellar stack equation of state” results from equating the water chem¬ 
ical potential in a lamellar stack of given composition (()-and, therefore, of given 
optimal interfacial area E((/))-submitted to the (hydrostatic) pressure p -|- H to 
the pure water chemical potential when pure water is submitted to (hydrostatic) 
pressure p. The resulting standard expression is 

n = - /.„ (6) 

either expressed as a n((/)) or a n(I') relation. 
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(a) Osmotic pressure 11 


(b) Lipid interfacial area S 


Figure 1: (a) Lamellar stack equation of state. Note the quasi-exponential 
decay of osmotic pressure with stacking period ^-continuous line. Decay length 
ca. 0.143 nm. (b) Bilayer equation of state. Dehydration results in a decrease in 
optimal interfacial area per lipid molecule. Redrawn from data for egg lecithin 
in Ref. [5] 

Examples of the two above-mentioned equations of state are given for illus¬ 
tration in Fig. [Tal and llbl-data extracted from Ref. [S]. The quasi-exponential 
decay of the osmotic pressure 11 with stacking period £ is also observed when 11 
is plotted against the thickness of the water channel £w = ^(1 — </>), though with 
a different decay length: A value ca. 0.256 nm is found in the latter case [5], 
instead of ca. 0.143 nm in Fig. Hal Quite significant also is the decrease in 
optimal interfacial area per lipid molecule E as dehydration proceeds, a feature 
already observed more than 50 years ago with surfactant-water or lipid-water 
systems [EHIZ]- 

In the quite common case where a phase coexistence between a lamellar 
structure with finite spacing £max and (almost) pure water when added in ex¬ 
cess is observed, the lamellar stack is said to be “bound”. This amounts to 
saying that, at some “dilution limit”, there is a finite value (p* for the lipid 
volume fraction where the osmotic pressure reaches zero: n((/)*) = 0. In the 
system studied in Ref. [S], for instance, the dilution limit is found at about 
^max = 6.25 nm. At contrast, (less common) “unbound” systems may be swollen 
seemingly indefinitely with water, and (()*—>■ 0 or, equivalently, £max —t -koo. 
The so-called “unbinding transition”, theoretically described in Ref. [TSl - ED] . 
separates “bound” systems, with somehow weak undulation interactions, from 
“unbound” systems where stronger undulation interactions overcome attractive, 
van der Waals forces. It has been experimentally evidenced in appropriately 
chosen lamellar systems, see for instance Ref. m for a recent report. 


2.2 Simple fluids 


As argued in more details in this Section, the above-described unbinding fea¬ 
tures of lamellar stacks upon addition of solvent has actually some similarities 
with the liquid-gas transition in pure compounds when hydrostatic pressure is 
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isothermally decreased: Molecular (translational) kinetic energy (3D configura¬ 
tional entropy) then plays the role of bilayer undulations in lamellar stacks (ID 
configurational entropy), with cohesive forces-viz. van der Waals interactions- 
acting similarly in both cases. Close to the triple point temperature, the liquid 
phase may be said to be “bound” in the sense that pressure may be decreased to 
(almost) zero while keeping a rather small specific volume for the system. Close 
to the critical point temperature, however, the liquid phase would be charac¬ 
terised as “unbound” because the specific volume may reach large values even 
though a significant pressure is still applied to the fluid. Fig. [5] illustrates this 
very common behaviour in the case of the fluid phase of ethane-data extracted 
from Ref. [HJ-where, for the sake of comparison with lamellar phase data, the 
physically-relevant parameter, namely the specific volume, has been expressed 
as an average distance i between neighbouring molecules in the liquid. 

A conspicuous feature of the experimental pressure equations of state dis¬ 
played in Fig.[2]may be noted at the high pressure end of the two curves, with a 
quasi-exponential behaviour p(£) oc exp(—f/A) that somehow resembles the one 
observed in Fig. [T^ At contrast with lamellar phase data, however, the char¬ 
acteristic decay length A is here very short. Values in the range 0.02-0.024 nm 
are found for A, the shorter value being appropriate for the lower temperature. 

Interpreting parameter A in terms of molecular dimensions seems difficult, 
but a simple interpretation for the observed behaviour is easily found in consid¬ 
ering the van der Waals equation of state for pure, fluid compounds. In terms 
of only two parameters, namely a second virial coefficient &2 describing the in- 
termolecular interactions, and an excluded volume vq accounting for the finite 
molecular dimensions, van der Waals equation of state is m 


2 Vo/V 


pvp _ 

ksT 2 Vq 



(7) 


V 


with V the specific volume of the van der Waals fluid, ks the Boltzmann constant 
and T the temperature. Dense fluids, i.e. systems with specific volumes only 
moderately larger than uq, may be obtained at moderate pressures if the second 
virial coefficient 62 is negative enough. With the integral expression for 62 in 
terms of the Mayer function [^, this occurs when fc^T, a measure of thermally- 
driven translational kinetic energy, is small enough compared to cohesive energy. 
Below a critical temperature where 62 decreases below the condensation 
threshold —27uo/4 in the van der Waals model, liquid phases may indeed be 
formed. An illustration is given in Fig. [3] for two values of the second virial 
coefficient, one close to, the other one farther below the condensation threshold. 
Note that, for plotting eq. ([7]) in Fig. [31 the specific volume range has been 
restricted: On the low density end, to avoid entering the liquid-gas coexistence 
region where eq. © no longer applies, but also on the high density end, to avoid 
a spurious logarithmic singularity implicit in the van der Waals description of 
the configurational entropy. 

The similarity between Fig. [2] and [3] is striking, in particular as regards the 
large mismatch between decay parameters, resulting from either experiment or 
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Figure 2: Pressure equations of state for liquid ethane extracted from Ref. m, 
expressed as a function of £ = j , at two temperatures T = 190 K 

(©) and T = 300 K (■). For pure ethane, triple point and, respectively, critical 
point temperatures and pressures are T* = 90.4 K and pt = 1.14 Pa, resp. 
Tc = 305.3 K and Pc = 4.93 MPa. The £ variable would be the nearest-neighbour 
distance if the mass density p of liquid ethane were obtained by placing molecules 
at the nodes of a face-centred cubic lattice. The molar mass of ethane is M = 
30.07 g/mol and Na is the Avogadro constant. Note the quasi-exponential decay 
of pressure for small intermolecular separations-decay lengths ca. 0.020 nm 
(T = 190 K) and 0.024 nm (T = 300 K) 
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Figure 3: Pressure equations of state in reduced units for the liquid phase of a 
van der Waals fluid. Second virial coefficients b 2 /vo = —10 (dashed line) and 
—7.92 (continuous line). From the common tangent construct, the liquid-gas 
binodal begins, respectively, at the specific volumes vl ~ 1.37uo and vl ~ 
1.67uo. Note the quasi-exponential decay of pressure with decay parameters, 
respectively, O.OSIvq and 0.143uo 
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model, and characteristic molecular dimensions. A compelling interpretation 
for the quasi-exponential pressure decay observed over a range of intermolec- 
ular separations is therefore the following: It is not related to any structural 
property of the system and should be attributed to the standard competition 
(with a temperature-dependent outcome, obviously) between cohesive energy, 
favoured in high density states because of attractive van der Waals interactions 
between molecules (62 term in eq. (O), and configurational entropy^ reduced in 
high density states owing to short-range steric repulsions (vq term in eq. Q). 
Besides, the precise shape of the interaction profile, e.g. a repulsive well 

combined with a —attractive component in the standard Lennard-Jones 
parametrisation of van der Waals interactions [23], is irrelevant here since the 
spatial integration of the Mayer function leading to 62 smoothes out such details. 

Taking one step further and considering the similarities between Fig. Hal with 
either Fig.|2|or Fig.|21 it now seems natural to relate the quasi-exponential decay 
in the lamellar stack equation of state to a similar competition between cohesive 
energy (arising in lamellar stacks from van der Waals interactions through the 
water channels between lipid bilayers) and configurational entropy, as repeatedly 
suggested by J. Israelachvili and H. Wennerstrom [M27j. This idea-in essence 
that “the swelling process is entropy driven” | 2 H]-has found a partial support 
in recent numerical simulations m- It is developed in Section |23l below. 

2.3 Lamellar stacks “a la van der Waals” 

Before presenting our model, a word of caution: Nothing prevents, in principle, 
some direct (repulsive) “hydration forces” playing a (marginal) role in the de¬ 
tailed characteristics of the lamellar stack equation of state. Possible physical 
mechanisms for such direct forces involve the solvent (water) binding properties 
with the bilayer surfaces m. or solvent structural properties [aiEa through 
a mechanism similar to wetting [83] in this latter case: With anti-symmetric 
boundary conditions at the bilayer-solvent interface, the order parameter pro¬ 
file of a (non scalar) property associated to solvent across the water channel, 
electric polarisation typically, would lead indeed to an effective repulsive inter¬ 
action in the lamellar stack. Interactions would nevertheless remain essentially 
attractive in the case of symmetric boundary conditions, as shown by Richetti 
et al. in their study of a structure-prone (pre-smectic) solvent confined between 
solid surfaces [33] ■ In any case, “hydration forces” if they exist in this restricted 
sense would merely affect the magnitude of the second virial coefficient in the 
van der Waals fluid analogy, increasing I 62 I (if attractive-more cohesion) or de¬ 
creasing 1621 (if repulsive-less cohesion), which would only slightly displace the 
location of the critical temperature in the phase diagram. 

It therefore appears necessary to reconsider carefully the interpretation of 
statements similar to 

Hydration repulsion universally acts between well-solvated surfaces in 

water and balances the van der Waals attraction in the nanometre range 


(quoted from Ref. [53]) found in references spanning more than 30 years, for in¬ 
stance Ref. 0111151. The underlying vision of a balance is commonly expressed 
in quantitative terms by adding pressures (forces per unit bilayer area), with 
as many contributions to the sum as identified interactions. In the simplest 
case where only hydration and van der Waals interactions are considered to be 
relevant, the total pressure is 0 


P{i)=Poexp[-{e-S)/X] 


H 

Ott 




{£-s)3 (e + s) 


2 


( 8 ) 


where the exponentially-decaying, repulsive term represents hydration forces 
and the combination of power laws describes, in the simplest possible way, at¬ 
tractive van der Waals forces between identical, planar and parallel objects of 
infinite lateral extension and thickness S [= 2v/Y, in eq. (|1|)], separated by a 
channel of width i — 6. In eq. (0), Pq > 0 is the amplitude of the hydration 
pressure, A its decay length and H > 0 the Hamaker coefficient. 

From Ref. 0, values appropriate for describing the experimental measure¬ 
ments in egg-lecithin-based lamellar stacks are: Pq ~ 7.05 x 10® Pa, A ~ 
0.256 nm and H Ri 6.0 x 10“^^ J. As already mentioned in Section 12.21 pa¬ 
rameter A-introduced here-differs from parameter A used for describing the 
lamellar stack equation of state n(£) displayed in Fig. [1^ This comes from the 
choice made in Ref. 0, the lamellar stack equation of state being represented 
by a different function W[£—5). It differs from n(£) because 5 actually depends 
on hydration. The “balance” between hydration and van der Waals forces is 
obtained at a maximum swelling fmax ~ 6.3 nm, where bilayer thickness is 
5 K, 3.5 nm. 

But there is more in thermodynamics than in mechanics or, in other terms, 
free energy is more than mere potential energy: A complete description of, e.g., 
equations of state or phase diagrams is not wholly encapsulated in force models 
similar to eq. (0) that ignore entropy, a criticism already explicitly formulated 
in the context of the unbinding transition by R. Lipowski and S. Leibler [I8j or 
S.T. Milner and D. Roux m- They showed that, even though the mechanical 
(but ad hoc) model, eq. 0), accounts extremely well for the measured osmotic 
pressure data-FiEf. llaH t utterly fails in explaining the detailed thermodynamic 
features of the unbinding transition. What is missing in the naive mechanical 
approach is, specifically, a description of the bilayer equation of state-Fisf. llbl -as 
well as a physically well-grounded interpretation of the “hydration forces”-two 
aspects recently emphasised in References [351[35] and [37], respectively. 

In the present approach, we propose to overcome these shortcomings by 
modelling the fundamental thermodynamic quantity in the problem, namely 
the excess free energy density /exc, with the following physical ingredients: i) 
Helfrich undulation (entropic) interactions between bilayers across the water 
channels, ii) van der Waals (and possibly other) direct interactions at the level 
of a second virial coefficient description as in Ref. [531 , and iii) bilayers described 
as being, individually, a kind of two-dimensional van der Waals fluid well below 
its own unbinding transition in the sense given to it in Section 12.21 It is worth 
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Table 1: Numerical values chosen for illustrating the properties of the thermo¬ 
dynamic model, eq. m 


mentioning again that, owing to the virial approach chosen in ii), the precise 
shape of the interaction potential, a repulsive exponentially-decaying term added 
to an attractive power-law contribution as would result from eq. ([S|) for instance, 
is irrelevant here, for the same reasons as given in Section 12.21 In quantitative 
terms, our heuristic approach amounts to writing 


^ ^ y*exc 


(ksT)^ 1 
m K (£ - (5)2 
-keTx^^e 

+ gln[ 

E Le(E-So) 



(9) 


where £ depends on E and (j) according to eq. o and S is again a convenient 
notation for 2w/E. In eq. ([51), e being the base of the natural logarithm, k 
is the bilayer bending modulus (controlling the amplitude of thermally-driven 
bilayer undulations 0), X is the second virial coefficient accounting for (inter¬ 
bilayer) direct interactions according to S.T. Milner and D. Roux m , Eo is 
the “excluded area” in the two-dimensional fluid of lipid molecules, &2 is the 
second virial coefficient accounting for intra-bilayer interactions between lipid 
molecules, “a la van der Waals”, and c a parameter with length units analogous 
to the thermal de Broglie wavelength in ideal (3D) classical gases |38) . 

For illustrating the essential role of ingredient iii) in the physical content of 
eq. Fig. m displays a set of curves for /exc as a function of the area E per 
lipid molecule for a few selected bilayer volume fractions 4>. The numerical val¬ 
ues for the parameters used in the computation of /exc are given in Table [TJ As 
suggested by the trends observed in Fig. 0] (and confirmed by a thorough numer¬ 
ical analysis), the “optimal” area where /exc reaches its minimum value depends 
only weakly on the bilayer volume fraction </> for dilute lamellar stacks, but much 
more strongly as the lamellar phase becomes highly dehydrated. This is easily 
interpreted, in qualitative terms, by a competition between stack entropy, as de¬ 
scribed by the Helfrich undulation term that strives for large I (corresponding 
to small E) values, and two-dimensional liquid bilayer entropy, as described by 
the excluded area Eq term, oppositely striving for large E values. The second 
virial coefficient 62 , negative enough to ensure a “bound” two-dimensional liquid 
state for the lipid fluid, comes into play when the stack entropy is optimised, 
i.e. for dilute enough lamellar stacks, and stabilises E to a value greater than, 
but still rather close to Eq (~ 1.4 x Eq for Table [1] parameter values). 

Even though it appears difficult to obtain the bilayer equation of state by 
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Figure 4: Excess free energy per unit volume of a bilayer stack /exc, according 
to the thermodynamic model eq. ([5]), drawn as a function of the area E per 
lipid molecule for a few selected values of the bilayer volume fraction: (p = 0.25- 
rightmost continuous line, (p = 0.80-dashed line, p = 0.90-dotted line, and 
p = 0.95-leftmost continuous line. Note that, for the sake of clarity, a in¬ 
dependent constant has been added to /exc to adjust the minimum of each curve 
to the (arbitrarily chosen) value 0.01 fc^r/nm^ (thin horizontal line). The shifts 
in free energy per unit volume are, along decreasing hydration, A si 0.46, 1.43, 
1.43 and 0.87 fc^T/nm^ 
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(a) Osmotic pressure H (b) Lipid interfacial area S 

Figure 5: (a) Osmotic pressure 11 of the lamellar stack as a function of the 
stacking period 1. The osmotic pressure reaches ca. 4 MPa when the lamellar 
phase is highly dehydrated. At the dilution limit fmax ~ 6.1 nm, 11 —>■ 0. 
There is a quasi-exponential decay of the osmotic pressure with a decay length 
K, 0.43 nm between these two limits, (b) Optimal area S per lipid molecule as 
a function of the bilayer volume fraction ij). Numerical results derived from the 
thermodynamic model eq. ([9|) with parameters given in Table [T] 


analytically solving eq. a numerical procedure for determining the optimal 
E((^) is easily implemented, considering the simple shapes of the /exc curves 
displayed in Fig. [H But the possibility of describing a “bound” lamellar stack, 
that is to say a system where the osmotic pressure 11 reaches 0 at a finite 
hydration (p* fSection 12.11) should be simultaneously taken into account. This 
is not a too stringent constraint since eq. (IHl) leads in the context of the present 
model eq. ([H]) to a rather simple expression for the osmotic pressure, namely 


Stt^ {ksTY 1 
64 K (£ — 




( 10 ) 


where the terms originating from the two-dimensional fluid of lipid molecules 
exactly cancel out. The resulting lamellar stack equation of state, as well as 
the bilayer equation of state are displayed in Fig. [Sal and [50 respectively. The 
broad features of both figures are rather similar to the experimental results 
displayed in Fig. [10 and Hbj for lamellar stacks, and (as far as the pressure 
equation of state is concerned) in Fig. [0 for a simple fluid. In addition, the 
unbinding transition remains well described within the framework of eq. (|9|). 
extending Ref. m- We have (numerically) explored the consequences on the 
dilution limit fmax of decreasing the Milner-Roux virial coefflcient y, keeping 
all other parameters fixed to the values given in Table [T] As shown in Fig. [0 
the prediction of our model is compatible with a scaling law ^max fx lx ~ Ac|~“) 
with critical unbinding occurring for Xc ~ 2.2 x 10“® nm“^ (or ksTxc ~ 9.2 Pa) 
and an unbinding exponent a ~ 0.9. The discrepancy between our result for the 
unbinding exponent, namely a ~ 0.9, and the theoretically-predicted values ip = 
1 in Ref. [5D] or ip = 1.00 ±0.03 in Ref. [T5] may come from the unsophisticated 
numerical methods used for obtaining the bilayer equation of state from eq. ©• 
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Figure 6: Dilution limit im&x as a function of \xlXc ~ 1|- The continuous line 
describes a power-law divergence of £max as the Milner-Roux virial coefficient 
approaches the critical value Xc, all other parameters in eq. being kept 
fixed to the values given in Table [TJ The power-law divergence, with ksTxc ~ 
9.2 Pa and a ~ 0.9, is the expected signature of a critical unbinding transition 
occurring in the lamellar stack as direct inter-bilayer interactions become less 
attractive [ID] 
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Figure 7: Small-angle x-ray scattering data for a series of lamellar stacks with 
bilayer composition x = 0.3 in osmotic equilibrium with polymer solutions of 
increasing concentrations: Cp = 0.068 g/mL (□), 0.100 (■), 0.141 (©), 0.162 
(•), 0.259 (A), 0.378 (a), and 0.472 (v). The stacking period i decreases with 
increasing osmotic pressure 


We leave a more detailed discussion to interested specialists. 

In Section [31 we further strengthen the relevance of our thermodynamic 
model, confronted with experimental results obtained for lecithin-Simulsol lamel¬ 
lar phases of varying bilayer compositions. 


3 Experimental results and discussion 

Small-angle x-ray scattering data are recorded, for a given bilayer composition 
X (described, from now on, in terms of Simulsol mass fraction, instead of mole 
fraction as in Section [T]), after immersing the dry lipid material in polymer 
solutions of various initial concentrations and waiting for water activity to equi¬ 
librate in both lamellar and polymer phases. Typical results are displayed in 
Fig. |7[ for a system with Simulsol mass fraction x = 0.3. The first and second 
order Bragg peaks, with ratio 1:2 as expected for a one-dimensional lamellar 
order, are always observed, a third order peak being also noticeable in the in¬ 
vestigated g-range for the more hydrated samples. The Bragg peaks are shifted 
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Figure 8: Osmotic pressure as a function of polymer concentration for three 
series of polymer solutions differing by the polymer degrees of polymerisation. 
PVPIO, A - PVP40, © - PVP360, □. The dotted line is a polynomial fit, 
appropriate for semi-dilute polymer solutions, with equation 11 = 4 • 10® x Cp x 
(1 + 6.25 X Cp) - pressure in Pa, concentration in g/mL 


towards larger values as osmotic pressure is increased because of the concomitant 
dehydration of the lamellar stack. 

From the first order Bragg peak location qo, the stacking period £ is directly 
obtained: £ = ‘1'k jq^. Since the dilution law, i.e. the £{^) relation, has been pre¬ 
viously experimentally determined |13j . the x-ray measurement indirectly gives 
the bilayer volume fraction (/. Because the amount of (initially dry) lipid in the 
whole system is known, mass conservation allows determining the equilibrium 
polymer concentration from its initial value in the known amount of polymer 
solution used for hydrating the lamellar stack. The lamellar stack equation of 
state then results from the auxiliary calibration curve (see Fig. [5]) giving the 
osmotic pressure FI of the polymer solution as a function of polymer concentra¬ 
tion. An illustration for the system with Simulsol mass fraction x = 0.3 is given 
in Fig. El with a quasi-exponential decay of the osmotic pressure observed in a 
range of stacking parameters £. 

The set of data shown in Fig.[T0]is obtained by repeating the same procedure 
and analyses for mass fractions spanning the pure lecithin (x = 0) to the almost 
pure Simulsol (x = 0.8) range. It was not possible to include the pure Simulsol 
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£ [nm] 

Figure 9: Lamellar stack equation of state n(£) for the system with bilayer 
composition x = 0.3. Continuous line from model discussed in Section 12.31 see 
text for details 
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Figure 10: Lamellar stack equations of state n(^) for systems with bilayer com¬ 
positions a; = 0 (□), X = 0.05 (■), x = 0.1 (©), x = 0.3 (•), x = 0.5 (A), x = 0.7 
(a) and a: = 0.8 (v). Continuous lines from model discussed in Section [531 see 
text for details 
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Figure 11: Experimental (symbols) and selected model (continuous lines) bilayer 
equations of state for lamellar stacks with bilayer compositions a: = 0 (□), 

X = 0.05 (■), X = 0.1 (©), X = 0.3 (•, —), X = 0.5 (A, —), x = 0.7 (A, — ) and 
X = 0.8 (V, —) 


system, and not even the x = 0.9 bilayer composition in the present study 
because at low osmotic pressures the water channel height reaches values so 
large that the polymer molecules can no longer be considered as excluded from 
the lamellar stack. Thus, increasing x apparently plays the role of approaching 
the unbinding transition here, as recently observed for a closely related system 
by a different technique M- 

The thermodynamic model discussed in Section 12.31 and, more specihcally, 
eq. (HI) have been used for describing the lamellar stack equations of state dis¬ 
played in Fig. El and [IHl as well as the bilayer equations of state displayed in 
Fig. [ni We have not attempted to implement a complete fitting procedure to 
finely tune values for all the parameters listed in Table (TJ but checked in a few 
numerical simulations the influence of changing the bilayer bending modulus k 
(bilayer stack entropy), as well as the Milner-Roux virial coefficient x (bilayer 
stack enthalpy, resulting from both attractive and repulsive contributions to 
inter-bilayer direct interactions) and excluded area Eg (two-dimensional fluid 
entropy). The continuous curves superimposed to data points in Fig. [9] and [TH] 
obviously give a fair description of the osmotic pressure data, with parameters 
given in Table [H The description, though acceptable, is less satisfactory as 
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X 

n/ksT 


So 

f 

■^max 



[nm-*^] 

[nm^] 

[nm] 

0.8 

0.55 

0.18 

0.370 

9.01 

0.7 

0.15 

0.43 

0.477 

8.86 

0.5 

0.047 

0.43 

0.635 

12.0 

0.3 

0.13 

0.43 

0.525 

8.53 

0.1 

0.33 

0.30 

0.495 

7.51 

0.05 

0.40 

0.36 

0.485 

6.95 

0.0 

0.35 

0.65 

0.505 

6.22 


Table 2: Numerical values chosen for describing experimental lamellar stack 
equations of state using eq. ([H])- Parameters v and & 2 /S 0 kept constant to values 
1.25 nm^ and —10, respectively. Parameter c = ^/T^. The last column gives 
the model-predicted dilution limits £max that result from the chosen parameters 


far as the bilayer equation of state is concerned, see Fig. [H] Note that, for a 
better readability, a (representative) subset only of the available model curves 
is displayed. Improvements in the two-dimensional liquid model for describing 
the intra-bilayer thermodynamic contribution to the lamellar stack excess free 
energy, eq. ©, might here be required. 

There is a rather remarkable trend in the evolution with the bilayer compo¬ 
sition X of the model parameters k, y and Sq) as well as in the model-predicted 
dilution limit £niax illustrated in Fig. I12al I12b[ I12cl and I12dl respectively. Two 
conspicuous compositions can be distinguished, one in the range x = 0.05-0.1 
close to pure lecithin, and when x is close to 0.5. In a previous structural study 
of the same lecithin-Simulsol system m, the value x = 0.5 was already exper¬ 
imentally evidenced as separating two regimes, even though the experimental 
dilution limit did not follow what is predicted here Fig. I12dl Still, the data for 
the area per (effective) lipid molecule S as a function of hydration displayed in 
Fig. [TT] clearly indicates that composition x = 0.5 experimentally plays a special 
role. The area per lipid molecule is significantly larger-and specially sensitive to 
hydration-than with all other bilayer compositions, an experimental feature in¬ 
dependent from any model. This could indicate that the composition threshold 
for the “brush-to-mushroom” conformational transition in the hydrophilic blocks 
of Simulsol molecules suggested in Ref. m and m is close to x = 0.5, and may 
give a clue for reaching a better agreement between experimental and predicted 
bilayer equations of state, at the expense of introducing a “brush-to-mushroom” 
variable into eq. ®. Besides, the experimental “back-and-forth” variation of 
S when x is close to zero observed in Fig. 1111 though not easily interpreted, is 
nicely compatible with the picture that emerges from the thermodynamic model 
of the lamellar stack, Fig. [T21 being simultaneously much less mo del-dependent: 
It should therefore correspond to a robust, but still unexplained, feature of the 
lecithin-dilute Simulsol mixture. 
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(a) Bending modulus k 



X 


(c) Excluded area Eq 



X 


(b) Milner-Roux second virial coefficient x 



(d) Dilution limit fmax 


Figure 12: Evolution as a function of bilayer composition x of the lamellar stack 
thermodynamic parameters. The continuous lines are guides for the eye. The 
bilayer configurational entropy is maximal close to x = 0.5-Fig. I12al Inter¬ 
bilayer attractive interactions are maximal and of comparable magnitude for 
X = 0 and x = 0.5-Fig. I12bl The excluded area of the two-dimensional bilayer 
fluid is maximal close to x = 0.5-Fig. I12cl The dilution limit of the lamellar 
stack is maximal close to x = 0.5-Fig. |l2dl 


4 Conclusion 

The interplay between soft confinement and steric effects has been studied in 
lamellar stacks of mixed lecithin-Simulsol bilayers, varying hydration and bi¬ 
layer composition, by coupling osmotic pressure with small-angle x-ray scat¬ 
tering measurements that give access to both “lamellar stack” and “bilayer” 
equations of states. 

With the help of a thermodynamic model of lamellar stacks that extends 
the approach of S.T. Milner and D. Roux by the explicit inclusion of the 
thermodynamic properties of a two-dimensional fluid representing the lipid bi¬ 
layers, we are able to quantitatively reproduce the osmotic pressure, and semi- 
quantitatively the lipid area data, keeping a meaningful description of the un¬ 
binding transition. The physical ingredients of the model are the Helfrich undu¬ 
lation interactions, controlled by a bending modulus parameter k, direct inter¬ 
bilayer interactions controlled by the Milner-Roux virial coefficient x, and a 
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description of the two-dimensional lipid fluid by the venerable van der Waals 
model designed more than a century ago for ordinary fluids, with an excluded 
area parameter Eq and a second virial coefficient & 2 - The celebrated, but elusive 
“hydration forces” play no role whatsoever in the model, while quasi-exponential 
osmotic pressure profiles remain, at least in a restricted range of water content 
in the lamellar stacks. 

Among the intriguing outcomes of applying the model to describing the 
available osmotic and small-angle scattering data, it appears that adding a co¬ 
surfactant (Simulsol) to the lipid bilayer not necessarily always increases the 
bilayer flexibility, even though k remains very sensitive to the amount of co¬ 
surfactant. This may be due to a kind of “brush-to-mushroom” transition taking 
place in the hydrophilic block of the Simulsol surfactant. The model prediction 
also does not quite fit with previous experimental results on the same system 
as far as the dilution limit of the lamellar stacks are concerned. 

The better knowledge of the lamellar stack free energy gained here should 
still be refined, for instance in view of understanding the detailed mechanisms 
of the lamellar-lamellar phase coexistence at low hydration observed in rather 
similar lecithin-Simulsol systems El- 
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A Appendix 

A.l Polymer solutions 

The polymer used in this study is polyvinylpyrrolidone (PVP), from Sigma- 
Aldrich, in three molecular weights (10, 40 and 360 kg/mol). The polymer 
solutions are purified using snakeskin dialysis tubes (from Thermo Scientific) 
of different pore dimensions. The dialysis compartments filled with polymer 
solution were immersed in a tank initially containing pure water, and the water 
content was changed a few times, until the total volume reached 300 times the 
volume of the solution. The whole procedure took about one week. Once pu¬ 
rified, the solutions were freeze-dried, which allowed re-dispersing the polymer 
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material in water at various desired concentrations. The complete polymer dis¬ 
solution could take a few days for the more concentrated solutions with PVP360. 

The osmotic pressure of each prepared solution was measured with the PZL- 
100 osmometer from PZL Tecnologia Company. Results are displayed in Fig. [SI 
The nearly-quadratic increase of osmotic pressure with polymer concentration, 
independently of the polymer molar mass, is characteristic of polymer solutions 
being in their semi-dilute regime. The n(cp) relation is empirically described 
by the polynomial law 11 = 4 • 10® x x (1 -|- 6.25 x Cp) - pressure in Pa, 
concentration in g/mL - which allows computing the osmotic pressure applied 
to lamellar stacks. 

A.2 Osmotically-stressed lamellar stacks 

The lipid membranes are prepared with soy lecithin (Avanti Polar Lipids) and 
a non-ionic commercial cosurfactant (Simulsol 2599 PHA, Seppic). Soy lecithin 
contains mainly dilinoleoylphosphatidylcholine (DLPC), with typically 35% of 
other (zwitterionic) lipids. Simulsol is a mixture of ethoxylated fatty acids 
(71% oleic and 11% palmitic acids being the main components), with an aver¬ 
age of 10 CH 2 -O-CH 2 groups in the hydrophilic block. Lecithin and Simulsol 
are co-solubilised in cyclohexane in desired proportions, which always leads to 
macroscopically homogeneous, transparent solutions. Therefore, a priori ho¬ 
mogeneous mixtures of controlled composition (labelled as a mass fraction x) 
are obtained by evaporating the solvent. Lamellar systems are then prepared by 
swelling a chosen mass of the dry amphiphilic mixture with thrice the amount of 
polymer solutions, tuning the polymer mass and initial concentration to scan as 
much as possible of the accessible domain of osmotic pressures once equilibrium 
is reached. The sample tubes were conserved at 4°C and cycles of centrifu¬ 
gation were done to accelerate homogenisation. Equilibrium is reached when 
water activity is the same in coexisting lamellar and polymer phases, which is 
easily deduced by observing the tubes characterised by a viscous and slightly 
turbid solution on top of a transparent polymer solution. Osmotic equilibrium 
also implies that some lecithin or Simulsol species are to be found in the poly¬ 
mer solution. The corresponding concentrations are, however, expected to be 
in the order of their respective critical micellar-aggregation-concentrations, viz. 
quite low. Neglecting the amount of “lost” material in computing the lamellar 
phase lipid fraction, as well as neglecting the contribution of lipid aggregates to 
osmotic pressure are therefore safe approximations. 

A.3 Small-angle x-ray scattering 

The small-angle scattering experiments were carried out at the Institute of 
Physics, University of Sao Paulo, Brazil, with the Xeuss instrument form Xenocs 
equipped with a Pilatus 300K detector (Dectris). Radiation produced by the 
microfocus Copper source is collected with a single-reflection multilayer op¬ 
tic producing a low-divergence, monochromatic beam with wavelength A = 
0.154 nm that is further collimated by a pair of scatterless slits-upstream slits 
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0.6x0.6 mm, downstream slits 0.5x0.5 mm. The sample-to-detector distance, 
calibrated with a Silver Behenate standard, is 0.77 m. Detector images span in 
practice a scattering wave vector range extending from 0.04 nm“^ to 3.5 nm“^. 
Samples are held in glass capillaries with a nominal diameter of 1.5 mm, and 
the scattering intensity is corrected for background by subtracting the properly 
normalised signal of a capillary filled with pure water. Exposure time varies 
from 15 up to 30 minutes, depending on the polymer solution concentration. 
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